arXiv: 1501.03395v2 [math.CO] 23 Apr 2015 


Calculation of the Number of all Pairs of Disjoint 
S-permutation Matrices 

Krasimir Yordzhev 


Faculty of Mathematics and Natural Sciences, South- West University , 2700 

Blagoevgrad, Bulgaria 
Email address: yordzhev@swu.bg 

Abstract: The concept of S-permutation matrix is considered. A general 
formula for counting all disjoint pairs of n 2 x n 2 S-permutation matrices as a 
function of the positive integer n is formulated and proven in this paper. To do 
that, the graph theory techniques have been used. It has been shown that to count 
the number of disjoint pairs of n 2 x n 2 S-permutation matrices, it is sufficient 
to obtain some numerical characteristics of all n x n bipartite graphs. 

Keyword: Binary matrix; S-permutation matrix; Sudoku matrix; Disjoint 
matrices; Bipartite graph 

MSC[2010]: 05B20; 05C50; 65F30 


1 Introduction and notation 


The text of the present paper is essentially extended and modified to the pre¬ 
vious articles of the author Eli, while we have corrected all inaccuracies and 
calculation errors. 

Let n be a positive integer. By [n] we denote the set [n] = {1,2,..., n}. 

A binary (or boolean , or (0,l)-matrix) is a matrix all of whose elements 
belong to the set f B = {0,1}. With we will denote the set of all n x n binary 
matrices. 

Let n be a positive integer and let A £ Q3„2 be a n 2 x n 2 binary matrix. With 
the help of n — 1 horizontal lines and n — 1 vertical lines A has been separated 
into n 2 of number non-intersecting n x n square sub-matrices Am, 1 < k,l < n, 
e.i. 

An A12 ■ ■ ■ Ai n 
A 21 A 22 ' • ' A2n 


A = 


( 1 ) 


A n 1 


A n 2 


A 


nn 


The sub-matrices Am, 1 < k,l < n will be called blocks. 
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A matrix A £ Q3 n 2 is called an S-permutation if in each row, each col¬ 
umn, and each block of A there is exactly one 1. Let the set of all n 2 x n 2 
S-permutation matrices be denoted by £ n 2 . 

As it is proved in [2] the cardinality of the set of all S-permutation matrices 
is equal to 

|S„ 2 | = (n!f\ (2) 

Two binary matrices A = (and B = ( bij ) € *8 m will be called 
disjoint if there are not elements with one and the same indices and such 
that CLij = = 1, i.e. if ciij = 1 then bij = 0 and if bij = 1 then = 0, 

1 < i,j < m. 

The concept of S-permutation matrix was introduced by Geir Dahl [2] in 
relation to the popular Sudoku puzzle. Sudoku is a very popular game. On the 
other hand, it is well known that Sudoku matrices are special cases of Latin 
squares in the class of gerechte designs T], 

Obviously a square n 2 x n 2 matrix M with elements of [n 2 ] = {1,2,..., n 2 } 
is a Sudoku matrix if and only if there are matrices A\, A 2 ,..., A n 2 £ £„ 2 , each 
two of them are disjoint and such that P can be given in the following way: 

M = 1 • A x + 2 ■ A 2 H-bn 2 - A n 2 (3) 

Some algorithms for obtaining random Sudoku matrices and their valuation 
are described in detail in [T]. 

In [2] Roberto Fontana offers an algorithm which randomly gets a family of 
n 2 x n 2 mutually disjoint S-permutation matrices, where n = 2,3. In n = 3 he 
ran the algorithm 1000 times and found 105 different families of nine mutually 
disjoint S-permutation matrices. Then using © he obtained 91-105 = 38 102 400 
Sudoku matrices. 

But it is known [4] that the total number of 9 x 9 Sudoku matrices is 
9! • 72 2 • 2 7 • 27 704 267 971 = 6 670 903 752 021 072 936 960 

Thus, in relation with Fontana’s algorithm, it looks useful to calculate the 
probability of two randomly generated S-permutation matrices to be disjoint. 
So the question of enumerating all disjoint pairs of S-permutation matrices nat¬ 
urally arises. This work is devoted to this task. 

Bipartite graph is the ordered triplet 

9 = {Rg U Cg , Eg), 

where R g and C g are non-empty sets such that R g D C g = 0, the elements of 
which will be called vertices. E g C R g x C g = {(r, c) | r £ R g , c £ C g } - the set 
of edges. Multiple edges are not allowed in our considerations. 

For more details on graph theory see M- 

Let n and k be integers such that 0 < k < n 2 . Let us denote with <3 n ,k the 
set of all bipartite graphs without multiple edges of the type g = ( R g ,C g ,E g }, 

n 2 

such that | R. g | = \C g \ = n and \E g \ = k. With 0 n = we denote the 

fc=0 
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set of all bipartite graphs without multiple edges of the type g = ( R g ,C g ,E g ), 
such that |i? 3 | = \C g \ = n and irrespective of the number of the edges. 

We will not take into consideration the nature of the vertices of the bipartite 
graphs, i.e. we will consider that for each of the examined bipartite graphs 
9 G ©n,fc> the vertices are respectively 

R g = {ri,r 2 , • ■ •, r n } 

and 

Cg — { C1 . C2, ■ • • , C n 

Let g' = (R g ' U C g >,E g >) and g" = (R g » L)C g »,E g ») be two bipartite graphs. 
We will say that g' and g" are isomorphic and we will write g' = g", if \R g ' \ = 
\R g " |, \C g /\ = \C g "\ and there exist bijections p : R g ' —> R g » and a : C g ' —> C g ", 
such that (r, c) € E g f •<==>■ (p(r),<r(c)) € E g ". 

Obviously the so entered relation = is an equivalence relation and with & n ,k 
we will denote the factor set (the set of the equivalence classes) 

©n,fc = / 

Let g = (R g ,C g , E g ) € © n> fc for some natural numbers n and k. Then with 
5 we will denote the set 


g = {he © n>fc | h = g} e ©„, fc 

of all isomorphic to g bipartite graphs, and with |p| their number, i.e. the 
cardinality of the set g. 

Let g = (R g ,C g , Eg) e & n ,k for some natural numbers n and k and let 
v e V g = R g U C g . With j(v) we will denote the set of all vertices from V g , 
adjacent to v, i.e. u e 7 (u) if and only if there exists an edge in E g which joins 
u and v. If v is an isolated vertex (i.e. there is no edge incident to v), then by 
definition j(v) = 0 and | 7 (f)| = 0. Obviously if v e R g , then j(v) C C gi and if 
v e C g , then 7 ( 11 ) C R g . 

Obviously 

1^)1 = 2n - 

vev g 

Let n and k be positive integers and let g e <8 n ,k■ We examine the ordered 
(n + l)-tuple 

W (g) = (V’o (g) ,ipi(g),...,ip n (g)), 

where tpi(g), i = 0 , 1 ,..., n is equal to the number of the vertices of g, incident to 
exactly i number of edges. It is easy to see that for each g e & n ,k the equalities 

n n 

^2ip z (g) = 2n, ^2iipi(g) = 2k. 

i=l i= 1 


have been executed. 
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For two bipartite graphs g and h, if g = h, then obviously 

W(sO = WO) = WO) = WO)- 

In the set ©„ we enter equivalence relation ” ~ ”, such that if A, B £ 23 n 
then A ~ B, if B is obtained from A after dislocation of some of the rows and/or 
columns of A. 

Let g = ( R g ,C g ,E g } £ ©„, where R g = {ri,r 2 , ? ..,r n } and C g = {ci,c 2 ,... ,c„}. 
Then we build the matrix A = [ay] G 23 n , such that ay = 1 if and only if 

(Vi ,Cj) £ Eg. 

Inversely, let A = [ay] £ 23„. We denote the i-th row of A with n, while the 
j-tli column of A with Cj. Then we build the bipartite graph g = (R g , C g ,E g ) £ 

©„, where R g = {ri, r 2 ,..., r„}, C g = {c\, c 2 ,..., c n } and there exists an edge 
from the vertex r,; to the vertex cq- if and only if ay = 1. 

Thus we showed the following obvious relation between the bipartite graphs 
and the binary matrices: 

Proposition 1 There exists one-to-one mapping 

T : 0 ra —> 

between the elements of & n and 03 such that if g,h £ & n , then 

9 — h <p(g) ~ <p(h). 

2 A representation of S-permutation matrices 

If Z\ Zn ... z n is a permutation of the elements of the set [n] = {1,2,... ,n} 
and we shortly denote p this permutation, then in this case we denote by p{i) 
the i-th element of this permutation, i.e. p(i) = Zi, i = 1 , 2 ,..., n. 

Let n„ denotes the set of all n x n matrices, constructed such that n £ U n 
if and only if the following three conditions are true: 

i) the elements of 7r are ordered pairs of numbers (i, j), where 1 < i,j < n; 

ii) if 

[<ai, 6i) (a 2 , bf) ••• (a n ,b n )] 

is the i -th row of 7r for any i £ [n] = {1,2 ,..., n}, then ai a 2 ... a n in this 
order is a permutation of the elements of the set [n] 

iii) if 

(ai,bi) 

(o 2 , bf) 
b n ) 

is the j-tli column of 7r for any j £ [n], then 6i,6 2 ,...,6 n in this order is a 
permutation of the elements of the set [n\. 
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From the definition it follows that every row and every column of any matrix 
of the set II n can be identified with permutation of elements of the set [n]. 
Conversely for every (2n)-tuple ((pi, p 2 , ■ • ■, p n ), (cri, cr 2 ,..., cr„)), where p* = 
Pi( 1) pi( 2) ... pi(n), <jj = Cj(l) <Tj( 2) ... <Jj(n ), 1 < i, j < n are permutations 
of elements of [n], then the matrix 


(pi(l),o-i(l)> (pi(2),cr 2 (l)) 

(P2(1),cti( 2)) (p 2 (2), ct 2 (2)) 

(Pn{ l),CTi(n)) (p n (2),a 2 (n)) 
is matrix of IL. Hence 


|n„| = (n!) 


2 n 


(pi(n),a n (l)} 

{p- 2 {n),a n { 2 )) 

(Pn(n), a n (n)) 


(4) 


Theorem 1 Let n be a positive integer, n > 2. Then there is one to one 
correspondence between the sets S „2 and H n . 


Proof. Let A G £ n 2 . Then A is constructed with the help of formula © 
and for every i,j G [n] in the block Aij there is only one 1 and let this 1 
has coordinates (a*, bj). For every i,j G [n] we obtain ordered pairs of numbers 
(a*, bj) corresponding to these coordinates. As in every row and every column of 
A there is only one 1, then the matrix [&ij} nxn , where aij = (a^, bj), 1 < i, j < n, 
which is obtained by the ordered pairs of numbers is matrix of n„, i.e. matrix 
for which the conditions i), ii) and iii) are true. 

Conversely, let [a»j] nxn G n n , where a = (a^bj), i,j G [n], ai,bj G [n]. 
Then for every i,j G [n] we construct binary n x n matrices Ajj with only one 
1 with coordinates (a*, bj). Then we obtain the matrix of type: ©■ According 
to the properties i), ii) and iii), it is obvious that the obtained matrix is S- 
permutation matrix. 

□ 

From Theorem [T|and formula (|U) it follows the proof of Proposition [3] in [2], 
i.e. formula ©. 

Definition 1 We say that matrices tt', tt" G H n , where tt' = \p'ij] n n > 7r " = 
[P"ij]nxn are disjoint, if p'ij 7 ^ p"^ for every pair of indices i,j G [n]. 

Proposition 2 The number of all pairs of disjoint matrices of £„2 is equal to 
the number of all pairs of disjoint matrices of A n . 

Proof. It is easy to see that with respect of the described in Theorem [l] one 
to one correspondence, every pair of disjoint matrices of £„2 will correspond to 
a pair of disjoint matrices of Tl n and conversely every pair of disjoint matrices 
of n n will correspond to a pair of disjoint matrices of £ n 2 . 

□ 


Definition 2 Let 7r',7r" G n n , n' = [//,,, 1 , tt" = and let the 

7 L r vJnxn 7 Ia UJnxn 

integers i,j G [n] are such thatpLj = p"\j. In this case we will say that p'^ and 
p"ij are component-wise equal elements. 
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Obviously two n n -matrices are disjoint if and only if they do not have 
component-wise equal elements. 

Example 1 We consider the following 113 -matrices: 



( 1 , 2 ) 

(3,1) 

(2,3) 

[p'ij\ = 

( 2 , 1 ) 

(3,3) 

( 1 , 2 ) 


(3,3) 

( 1 , 2 ) 

( 2 , 1 ) 


' (1,3) 

(3,2) 

( 2 , 1 ) 

m = 

(3,1) 

( 1 , 1 ) 

( 2 , 2 ) 


. (3,2) 

(1,3) 

(2,3) 


' ( 1 , 2 ) 

(3,3) 

( 2 , 1 ) 

[p'ij] = 

( 2 , 1 ) 

(3,2) 

( 1 , 2 ) 


(3,3) 

( 1 , 1 ) 

(2,3) 


Matrices tv' and tv" are disjoint, because they do not have component-wise 
equal elements. 

Matrices tv" and tv'" are not disjoint, because they have two component-wise 
equal elements: p' l3 = p" 3 = (2,1) and P 33 = p 33 = (2, 3). 

Matrices tv' and tv'" are not disjoint, because they have four component-wise 
equal elements: p' n = p'{[ = ( 1 , 2 ), p' 21 = p'f[ = ( 2 , 1 ), p'. 23 = p'J,' 3 = ( 1 , 2 ) and 

P31 = P31 = ( 3 , 3 )- 

3 A formula for counting disjoint pairs of n 2 x n 2 
S-permutation matrices 

Lemma 1 Let tv £ II„. Then the number q{n,k) of all matrices tv' £ II n , 
having at least k, k = 0, 1 ,..., n 2 component-wise equal elements to the matrix 
tv is equal to 

q{n, k) = Iffl 

Proof. Let tv = [Pij] nxn , tt' = [Pjj] £ n n and let tv and tv' have exactly 
k component-wise equal elements. Then we uniquely obtain the binary n x n 
matrix A = [dij] nxn i such that = 1 if and only if = pb, i,j £ [n]. Let 
graph g £ & n> k be such that g = p~ l {A), where <p : 03„ —> ©„ is the one-to-one 
mapping which gave us the grounds to formulate Proposition [L] This graph g 
identically corresponds to the ordered pair of matrices {n, tv') £ II x II. 

Inversely, let g = ( R g ,C g ,E g ) £ ©„ jfe , V g = R g U C g , A = [a i:j ] nxn = p(g) 
and let tv = \Pij\ nxn be an arbitrary matrix from II ra . We search for the number 
h{ tv, A) of all matrices tv' = [pT ] nxn £ II n , such that pL = Pij, if = 1. (It is 
assumed that there exist s, t £ [n] such that a s t = 0 and p' st = p st .) Let the i-th 
row, i = 1, 2,..., n of n correspond to the permutation pi of the elements of [n] 


/ n—2 


n k™ - on 


V’i(s) 


(5) 
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and let the ?’-th row of the matrix A correspond to the vertex r* £ R g of graph 
. Then there exist (n — ( 7 ( 77 ) |)! permutations p[ of the elements of [n], such 
that if an = 1, then pi(t) = p'(f), t £ [n]. Likewise we also prove the respective 
statement for the columns of it. Therefore 

h(n,A) = (n- 

vev g 

From everything said so far it follows that for each n £ II n there exist 


q(n,k) = Y n(»~ 

g£<8n,k \veVg 

matrices from II„, which have at least k elements that are component-wise equal 
to the respective elements of n. 

But obviously, if h £ is such a bipartite graph that h = g, then 


n(»- = n ( n - i^wd 1 ) = (n ^ 

iVEVg ) \u£V h ) \i—0 

whence follows the equality 


i’ii.g) 


q{n,k)= Y Iff! l( IT ( n ~l7(T)|)! = Y M II [(»-<)«] 


i’iia) 


gG&n,k \y£RgUCg 


gee n ,k 0 


And since (n — n)\ = 0! = 1 and [n — (n — 1)]! = 1! = 1, then we finally 
obtain formula 0 . 

□ 


Theorem 2 Let A £ E n 2 . Then the number of all matrices B £ E „2 which 
are disjoint with A is equal to 


tin 


( n —2 


(n!) 2n + ^(-l) fe Y Iffl II [("-W 


^i(g) 


k =1 




vi=0 


( 6 ) 


Proof. Let n > 2 be an integer. Then applying theorem [l] lemma [l] and the 
principle of inclusion and exclusion we obtain that the number of all matrices 
B £ E „2 which are disjoint with A is equal to 


fc=l 

where the function q(n,k ) is calculated with the help of formula (0, while |n„| 
with the help of formula 0. Thus we obtain the proof to formula 0. 

□ 
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Corollary 1 The cardinality rj n of the set of all disjoint non-ordered pairs of 
n 2 x n 2 S-permutation matrices is equal to 


Vn ~ ^ sr* 


(7) 


where is described using formula OJ 


Proof, follows directly from formula @ and having in mind that the ’’dis¬ 
joint” relation is symmetric and antireflexive. 

□ 

Corollary 2 The probability p{n) of two randomly generated n 2 x n 2 S-permutation 
matrices to be disjoint is equal to 


P(n) = 




(n\) 2n — 1 ’ 


where is described using formula [5J 

Proof. Applying Corollary [T] and formula ©, we obtain: 


(?r!) 


\\2n 


p{n) = 


<n 




i\) 2n ((n!) 


- 1 


( 8 ) 


(n\) 2n — 1 


□ 


4 Calculation of the Number of the Disjoint Pairs 
of S-permutation matrices when n = 2 and 

n = 3 

4.1 Consider n — 2 

When n = 2, ©2 is composed of seven equivalence classes also including the 
graph without edges (k = 0), which does not participate in our calculations. 
When k = 1,2, 3,4, we have depicted one representative g i, 52 , 53 , 54 , 55 and 
ge from each equivalence class respectively on figures [U [H [3] and [I] 

It is not difficult to notice that 
lffil=4, (V’Xffi) = <2,2,0), k= 1 
|P 2 |=2, (il>)(g 2 ) = { 0,4,0), k = 2 

|5 3 I = 11741 = 2, (if){g 3 ) = (ip)(g 4 ) = {1,2,1), k = 2 
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Figure 5: n = 3, k = 1 

l&l =4, mg 5 ) = ( 0,2,2), fc = 3 

i^ei = W^9&) = (0) 0j 4 ); k = 4 

Let j4 6 S4. Then the number £2 of all matrices B £ Y >4 which are disjoint 
with A is equal to 


6 = (2!) 4 + £(-l) fc ( £ \g\.2^)] = 

\gA& n , k 


k =1 


= 16 - 4 • 2 2 + (2 • 2° + 2 • 2 1 + 2 • 2 1 ) - 4 • 2° + 1 • 2° = 7 


2 4 

m = = 56 


P( 2) = 



7 

15 


4.2 Consider n = 3 

On figures from [5] to [13] one representative from each equivalence class of the 
factor sets ©3,/c, k = 1,2,..., 9 has been depicted, and we have numbered these 
graphs from 7 to 41. 

When n = 3 formula ([HI) has the form of 

41 

£3 = 6 6 + 15*16^2^) (-1)*®*) , 

i=7 

where with K,('g i ) we have denoted the number of the edges of the graphs from 
the equivalence class g,. 

Below we enumerate the examined numerical characteristics of the respec¬ 
tive equivalence classes of bipartite graphs. Their calculation is trivial. 

ls 7 l=(J) =9, W(ff 7 ) = < 4,2,0,0), k= 1 
l&l = (J) l&l = 18, <V'}(3s) = <2,4,0,0), fc = 2 
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Figure 6: n = 3, k = 2 
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Figure 7: n = 3, k = 3 
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Figure 8: n = 3, k = 4 
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Figure 9: n = 3, k = 5 




Figure 10: n — 3, k = 6 


13 













1091 = 10101 = (i) (9) = 9 ’ = (V’)(0io) = (3,2,1,0), k = 2 

I011I = |<S 3 | = 3! = 6, (V>)(0n) = (O,6,O,O), k = 3 

10121 = 10131 = (^j |03l = 18 > W(012) = W(013) = (M. 1 * 0 ). ^ = 3 

I014I = (j) |0 5 | = 36, W( 5 14 ) = (2,2,2,0), k = 3 

I015I = l0iel = ^ =3, W(0i5) = W(0ie) = ( 2 ,3,0,1), fc = 3 

I017I = (l) = 9 > (V J )(0i?) = (2,0,4,0), k = 4 

l0!8l=(i) ( 2 ) =9. W( 018 ) = (O,4,2,O), fc = 4 

I 01 9 1 = ( 2 ) I 05 I = 36, W(gw) = (0,4, 2 , 0 ), A; = 4 

l 02 ol = I 021 I = I 02 I = 13, ('0)(02o) = W( 02 i) = (1,2,3,0), fc = 4 

I 022 I = I 023 I = = 18, (^>( 022 ) = W(023) = (1,3,1,1), fc = 4 

10241 = = 9, (V’)(024) = (9,4,0,2), k = 5 

10251=^ =9, W ( 025 ) = (0,2,4,0), fc = 5 

10261 = ^ 2 ^) =36, (V , )(026) = (0,2,4,0), k = 5 

I 027 I = 10281 = ^ I 02 I = 18, W»)( 027 ) = (V’)( 02 s) = (0,3,2,1), fc = 5 

10291 = l03ol = 3! = 18, (^)(g 2 g) = (^)(03o) = (1,1,3,1), k = 5 

I 031 I = I 011 I = 6, (V’)(0 31 ) = (0,0,6,0), k = 6 

1032 I = 10331 = I 012 I = 18. W(032) = W(033) = (0, 1,4,1), k = 6 

10341 = I 014 I = 36, {ip)(g 34 ) = (0,2,2,2), fc = 6 

10351 = 10361 = 10151 = 3, W(035) = (V’)(036) = (1,0,3,2), fc = 6 

10371 = |0 8 I = 18, W(037) = (O,O,4,2), fc = 7 

10381 = l039l 2 = I09I = 9, (^)(038) = W(0 39 ) = (0,1,2,3), A; = 7 

|0 4O I= Q =9, (V’)(04o) = (0,0,2,4), fc = 8 

|04il = l> W(04i) = (0,0,0,6), fc = 9 

In order to calculate £3 we used the next computer programme written in 
programming language C++. 

int mainO { 
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double xi_3 =pow(6,6); 
int g[] [4] = { 

9,4,2,1, 

18.2.4.2, 

9.3.2.2, 

9,3,2,2, 

6,0,6,3, 

18,1,4,3, 

18.1.4.3, 

36.2.2.3, 

3.2.3.3, 

3,2,3,3, 

9,2,0,4, 

9,0,4,4, 

36,0,4,4, 

18,1,2,4, 

18.1.2.4, 

18.1.3.4, 

18.1.3.4, 

9,0,4,5, 

9,0,2,5, 

36,0,2,5, 

18,0,3,5, 

18,0,3,5, 

18.1.1.5, 

18,1,1,5, 

6 , 0 , 0 , 6 , 

18,0,1,6, 

18,0,1,6, 

36,0,2,6, 

3,1,0,6, 

3,1,0,6, 

18,0,0,7, 

9,0,1,7, 

9,0,1,7, 

9,0,0,8, 

1,0,0,9 

>; 

for (int i=0; i<35; i++) 

xi_3 += g [i] [0] * pow(6,g[i] [1] ) * pow(2,g[i] [2]) * (g[i][3]°/„2 ? -1 : 1) 
double p3; 
long eta_3; 

eta_3 = (long) (pow(6,6) /2) * xi_3; 
p3 =xi_3 / (pow(6,6)-1); 
cout<<"xi_3 = "<<xi_3<<endl; 
cout<<"eta_3 = "<<eta_3<<endl; 


16 



cout<<"p(3) = "<<p3<<endl; 
return 0; 

> 


After its work we obtain that 


& = 17 972 


Then according to formula 0, the number of all disjoint nonordered pairs 
of 9 x 9 S-permutation matrices is equal to 

(3!) 6 

r/ 3 = ^-£3 = 419 250 816. 

The probability p(3) for two randomly obtained 9x9 S-permutation matrices 
to be disjoint, according to formula 0 is equal to 

p( 3) =-^-= 0.385211 

V ' (3!) 6 — 1 
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